Modeling HER2 Effects on Cell Behavior from Mass Spectrometry Phosphotyrosine Data

Cellular behavior in response to stimulatory cues is governed by information encoded within a complex intracellular signaling network. An understanding of how phenotype is determined requires the distributed characterization of signaling processes (e.g., phosphorylation states and kinase activities) in parallel with measures of resulting cell function. We previously applied quantitative mass spectrometry methods to characterize the dynamics of tyrosine phosphorylation in human mammary epithelial cells with varying human epidermal growth factor receptor 2 (HER2) expression levels after treatment with epidermal growth factor (EGF) or heregulin (HRG). We sought to identify potential mechanisms by which changes in tyrosine phosphorylation govern changes in cell migration or proliferation, two behaviors that we measured in the same cell system. Here, we describe the use of a computational linear mapping technique, partial least squares regression (PLSR), to detail and characterize signaling mechanisms responsible for HER2-mediated effects on migration and proliferation. PLSR model analysis via principal component inner products identified phosphotyrosine signals most strongly associated with control of migration and proliferation, as HER2 expression or ligand treatment were individually varied. Inspection of these signals revealed both previously identified and novel pathways that correlate with cell behavior. Furthermore, we isolated elements of the signaling network that differentially give rise to migration and proliferation. Finally, model analysis identified nine especially informative phosphorylation sites on six proteins that recapitulated the predictive capability of the full model. A model based on these nine sites and trained solely on data from a low HER2-expressing cell line a priori predicted migration and proliferation in a HER2-overexpressing cell line. We identify the nine signals as a “network gauge,” meaning that when interrogated together and integrated according to the quantitative rules of the model, these signals capture information content in the network sufficiently to predict cell migration and proliferation under diverse ligand treatments and receptor expression levels. Examination of the network gauge in the context of previous literature indicates that endocytosis and activation of phosphoinositide 3-kinase (PI3K)-mediated pathways together represent particularly strong loci for the integration of the multiple pathways mediating HER2′s control of mammary epithelial cell proliferation and migration. Thus, a PLSR modeling approach reveals critical signaling processes regulating HER2-mediated cell behavior

Receptor Tyrosine Kinases Fall into Distinct Classes Based on Their Inferred Signaling Networks

Although many anticancer drugs that target receptor tyrosine kinases (RTKs) provide clinical benefit, their long-term use is limited by resistance that is often attributed to increased abundance or activation of another RTK that compensates for the inhibited receptor. To uncover common and unique features in the signaling networks of RTKs, we measured time-dependent signaling in six isogenic cell lines, each expressing a different RTK as downstream proteins were systematically perturbed by RNA interference. Network models inferred from the data revealed a conserved set of signaling pathways and RTK-specific features that grouped the RTKs into three distinct classes: (i) an EGFR/FGFR1/c-Met class constituting epidermal growth factor receptor, fibroblast growth factor receptor 1, and the hepatocyte growth factor receptor c-Met; (ii) an IGF-1R/NTRK2 class constituting insulin-like growth factor 1 receptor and neurotrophic tyrosine receptor kinase 2; and (iii) a PDGFRβ class constituting platelet-derived growth factor receptor β. Analysis of cancer cell line data showed that many RTKs of the same class were coexpressed and that increased abundance of an RTK or its cognate ligand frequently correlated with resistance to a drug targeting another RTK of the same class. In contrast, abundance of an RTK or ligand of one class generally did not affect sensitivity to a drug targeting an RTK of a different class. Thus, classifying RTKs by their inferred networks and then therapeutically targeting multiple receptors within a class may delay or prevent the onset of resistance.W. M. Keck FoundationNational Institutes of Health (U.S.) (R21 CA126720)National Institutes of Health (U.S.) (P50 GM068762)National Institutes of Health (U.S.) (RC1 HG005354)National Institutes of Health (U.S.) (U54-CA112967)National Institutes of Health (U.S.) (R01-CA096504)Alfred and Isabel Bader (Fellowship)Jacques-Emile Dubois (fellowship

Maximum Entropy Reconstructions of Dynamic Signaling Networks from Quantitative Proteomics Data

Advances in mass spectrometry among other technologies have allowed for quantitative, reproducible, proteome-wide measurements of levels of phosphorylation as signals propagate through complex networks in response to external stimuli under different conditions. However, computational approaches to infer elements of the signaling network strictly from the quantitative aspects of proteomics data are not well established. We considered a method using the principle of maximum entropy to infer a network of interacting phosphotyrosine sites from pairwise correlations in a mass spectrometry data set and derive a phosphorylation-dependent interaction network solely from quantitative proteomics data. We first investigated the applicability of this approach by using a simulation of a model biochemical signaling network whose dynamics are governed by a large set of coupled differential equations. We found that in a simulated signaling system, the method detects interactions with significant accuracy. We then analyzed a growth factor mediated signaling network in a human mammary epithelial cell line that we inferred from mass spectrometry data and observe a biologically interpretable, small-world structure of signaling nodes, as well as a catalog of predictions regarding the interactions among previously uncharacterized phosphotyrosine sites. For example, the calculation places a recently identified tumor suppressor pathway through ARHGEF7 and Scribble, in the context of growth factor signaling. Our findings suggest that maximum entropy derived network models are an important tool for interpreting quantitative proteomics data

Linear combinations of docking affinities explain quantitative differences in RTK signaling

Receptor tyrosine kinases (RTKs) process extracellular cues by activating a broad array of signaling proteins. Paradoxically, they often use the same proteins to elicit diverse and even opposing phenotypic responses. Binary, ‘on–off' wiring diagrams are therefore inadequate to explain their differences. Here, we show that when six diverse RTKs are placed in the same cellular background, they activate many of the same proteins, but to different quantitative degrees. Additionally, we find that the relative phosphorylation levels of upstream signaling proteins can be accurately predicted using linear models that rely on combinations of receptor-docking affinities and that the docking sites for phosphoinositide 3-kinase (PI3K) and Shc1 provide much of the predictive information. In contrast, we find that the phosphorylation levels of downstream proteins cannot be predicted using linear models. Taken together, these results show that information processing by RTKs can be segmented into discrete upstream and downstream steps, suggesting that the challenging task of constructing mathematical models of RTK signaling can be parsed into separate and more manageable layers

Effects of HER2 overexpression on cell signaling networks governing proliferation and migration

Although human epidermal growth factor receptor 2 (HER2) overexpression is implicated in tumor progression for a variety of cancer types, how it dysregulates signaling networks governing cell behavioral functions is poorly understood. To address this problem, we use quantitative mass spectrometry to analyze dynamic effects of HER2 overexpression on phosphotyrosine signaling in human mammary epithelial cells stimulated by epidermal growth factor (EGF) or heregulin (HRG). Data generated from this analysis reveal that EGF stimulation of HER2-overexpressing cells activates multiple signaling pathways to stimulate migration, whereas HRG stimulation of these cells results in amplification of a specific subset of the migration signaling network. Self-organizing map analysis of the phosphoproteomic data set permitted elucidation of network modules differentially regulated in HER2-overexpressing cells in comparison with parental cells for EGF and HRG treatment. Partial least-squares regression analysis of the same data set identified quantitative combinations of signals within the networks that strongly correlate with cell proliferation and migration measured under the same battery of conditions. Combining these modeling approaches enabled association of epidermal growth factor receptor family dimerization to activation of specific phosphorylation sites, which appear to most critically regulate proliferation and/or migration

2D protrusion but not motility predicts growth factor–induced cancer cell migration in 3D collagen

Growth factor–induced migration is a critical step in the dissemination and metastasis of solid tumors. Although differences in properties characterizing cell migration on two-dimensional (2D) substrata versus within three-dimensional (3D) matrices have been noted for particular growth factor stimuli, the 2D approach remains in more common use as an efficient surrogate, especially for high-throughput experiments. We therefore were motivated to investigate which migration properties measured in various 2D assays might be reflective of 3D migratory behavioral responses. We used human triple-negative breast cancer lines stimulated by a panel of receptor tyrosine kinase ligands relevant to mammary carcinoma progression. Whereas 2D migration properties did not correlate well with 3D behavior across multiple growth factors, we found that increased membrane protrusion elicited by growth factor stimulation did relate robustly to enhanced 3D migration properties of the MDA-MB-231 and MDA-MB-157 lines. Interestingly, we observed this to be a more reliable relationship than cognate receptor expression or activation levels across these and two additional mammary tumor lines.National Institutes of Health (U.S.) (Grant no. R01- GM081336)National Science Foundation (U.S.). Graduate Research FellowshipUnited States. Dept. of Defense. Congressionally Directed Medical Research Programs. Breast Cancer Research Program (Grant no. W81XWH-11-1-0088)United States. Dept. of Defense. Congressionally Directed Medical Research Programs. Breast Cancer Research Program (Grant no. W81XWH-10-1-0040

Current challenges in software solutions for mass spectrometry-based quantitative proteomics

This work was in part supported by the PRIME-XS project, grant agreement number 262067, funded by the European Union seventh Framework Programme; The Netherlands Proteomics Centre, embedded in The Netherlands Genomics Initiative; The Netherlands Bioinformatics Centre; and the Centre for Biomedical Genetics (to S.C., B.B. and A.J.R.H); by NIH grants NCRR RR001614 and RR019934 (to the UCSF Mass Spectrometry Facility, director: A.L. Burlingame, P.B.); and by grants from the MRC, CR-UK, BBSRC and Barts and the London Charity (to P.C.

MCAM: Multiple Clustering Analysis Methodology for Deriving Hypotheses and Insights from High-Throughput Proteomic Datasets

Advances in proteomic technologies continue to substantially accelerate capability for generating experimental data on protein levels, states, and activities in biological samples. For example, studies on receptor tyrosine kinase signaling networks can now capture the phosphorylation state of hundreds to thousands of proteins across multiple conditions. However, little is known about the function of many of these protein modifications, or the enzymes responsible for modifying them. To address this challenge, we have developed an approach that enhances the power of clustering techniques to infer functional and regulatory meaning of protein states in cell signaling networks. We have created a new computational framework for applying clustering to biological data in order to overcome the typical dependence on specific a priori assumptions and expert knowledge concerning the technical aspects of clustering. Multiple clustering analysis methodology (‘MCAM’) employs an array of diverse data transformations, distance metrics, set sizes, and clustering algorithms, in a combinatorial fashion, to create a suite of clustering sets. These sets are then evaluated based on their ability to produce biological insights through statistical enrichment of metadata relating to knowledge concerning protein functions, kinase substrates, and sequence motifs. We applied MCAM to a set of dynamic phosphorylation measurements of the ERRB network to explore the relationships between algorithmic parameters and the biological meaning that could be inferred and report on interesting biological predictions. Further, we applied MCAM to multiple phosphoproteomic datasets for the ERBB network, which allowed us to compare independent and incomplete overlapping measurements of phosphorylation sites in the network. We report specific and global differences of the ERBB network stimulated with different ligands and with changes in HER2 expression. Overall, we offer MCAM as a broadly-applicable approach for analysis of proteomic data which may help increase the current understanding of molecular networks in a variety of biological problems.National Institutes of Health (U.S.) (NIH-U54-CA112967 )National Institutes of Health (U.S.) (NIH-R01-CA096504

Phosphoproteomics-Based Modeling Defines the Regulatory Mechanism Underlying Aberrant EGFR Signaling

BACKGROUND: Mutation of the epidermal growth factor receptor (EGFR) results in a discordant cell signaling, leading to the development of various diseases. However, the mechanism underlying the alteration of downstream signaling due to such mutation has not yet been completely understood at the system level. Here, we report a phosphoproteomics-based methodology for characterizing the regulatory mechanism underlying aberrant EGFR signaling using computational network modeling. METHODOLOGY/PRINCIPAL FINDINGS: Our phosphoproteomic analysis of the mutation at tyrosine 992 (Y992), one of the multifunctional docking sites of EGFR, revealed network-wide effects of the mutation on EGF signaling in a time-resolved manner. Computational modeling based on the temporal activation profiles enabled us to not only rediscover already-known protein interactions with Y992 and internalization property of mutated EGFR but also further gain model-driven insights into the effect of cellular content and the regulation of EGFR degradation. Our kinetic model also suggested critical reactions facilitating the reconstruction of the diverse effects of the mutation on phosphoproteome dynamics. CONCLUSIONS/SIGNIFICANCE: Our integrative approach provided a mechanistic description of the disorders of mutated EGFR signaling networks, which could facilitate the development of a systematic strategy toward controlling disease-related cell signaling